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Представлены результаты математического и экспериментального моделирования процессов распростране- 
ния угарного газа и избытков теплоты в вВОЗДУшной среде гаражного помещения. На основе теоретико- 
экспериментального метода получены адекватные экспериментальным данным поля подвижности воздуха, 
концентрации оксидов углерода и температуры внутри исследуемого помещения. 
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Введение. Снижение эмиссии вредных веществ (ВВ) в рабочую зону операторов участков обкат- 
ки машиностроительных предприятий является актуальной задачей [1]. С другой стороны, выбро- 
сы при горячей обкатке машин «тёплые», что приводит к образованию в загазованной рабочей 
зоне областей с неудовлетворительными параметрами микроклимата. В результате формируется 
целый комплекс опасных и вредных производственных факторов (ОВПФ), таких как загазован- 
ность, повышенная подвижность воздуха и температура. Каждый из ОВПФ необходимо учитывать, 
так как они оказывают определённое влияние друг на друга, что в конечном итоге сказывается на 
точности моделирования. 

Постановка задачи. Математическое моделирование процессов распространения вредных 
веществ в воздушной среде вентилируемых помещений связано с рядом трудностей. 

Во-первых, результаты физического моделирования исследуемых процессов позволяют 
уточнить как сами модельные уравнения, так и их начальные и граничные условия. 

Во-вторых, доминирующее влияние на формирование полей загазованности, подвижности 
воздуха и температуры оказывают коэффициенты турбулентной диффузии (обмена) и 
теплопроводности газовой смеси в помещении, входящие в систему дифференциальных 
уравнений. В системе используются дифференциальные уравнения, описывающие процессы 
конвективно-диффузионного тепломассопереноса, введенные в работах М.Я. Берлянда [2], Г.И. 
Марчука [3], С.В. Патанкара [4]. 

Поэтому вполне обоснованной выглядит поставленная в работе задача — исследовать 
поведение математической модели, описывающей распространение угарного газа и избытков 
теплоты в воздушной среде заранее выбранного модельного объекта с простой геометрической 
формой и известным внутри стационарным источником загрязнений. 

В качестве такого объекта — газовой камеры — был выбран типовой металлический гараж. 
Источником загрязнений являлся автомобиль ВАЗ 2104, отработанные газы от которого были 
отведены в центр гаражного помещения. 
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Математическая модель тепломассопереноса вредных веществ и избытков явной теп- 
лоты. 
Система уравнений в общем виде 
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С, — концентрация 1-го ВВ в помещении, мг/м?; #-— текущее время, с; и, и, ии- компоненты 


вектора скорости воздушной среды, м/с; и — скорость оседания/подъёма примеси/газов; 
А’ =(А,,А,,А,) -— коэффициенты турбулентной диффузии в помещении, м’/с; 


Р=(р,,О,,О,) - коэффициенты диффузии 1-го газа, м/с; О/(®) — интенсивность выбро- 
са /-го вредного вещества от /-го источника в помещении, кг/с; б(’—и) — дельта-функция 
Дирака; г, = (х,,у,,2,) — координаты источника выброса, м; р -— плотность газовой смеси, 


кг/м; С ‚ - удельная изобарная теплоемкость, Дж/(кг ` К); ^, — коэффициент теплопровод- 
ности Г-го газа, Вт/(м'К); ^, — коэффициент турбулентной теплопроводности, Вт/(м'К); 


^› — коэффициент радиационной теплопроводности, Вт/(м-К); О/ - интенсивность внут- 


ренних источников тепла, Вт; Р, - диффузионное число Прандтля, Р, — турбулентное 


число Прандтля. 
Граничные условия для уравнения (1): 
- на боковых стенках: 
56 0; (5) 
ди 
в последнем равенстве указывается, что стенки помещения, в котором проводятся исследования, 
являются непроницаемыми для концентрации загрязняющего вещества; 
- на входной и выходной границе: 
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где «| м ‚ это означает, что поток примеси зависит от её концентрации, скорости о, и 


направления воздушного потока, индуцированного работой вентиляции. 
Граничные условия для уравнения (2): 
- на твердых стенках 
ОР _ 
дп 
где л — единичный вектор внешней нормали. Это означает, что стенки исследуемого помещения 
являются непроницаемыми для воздушных потоков; 
- на выходной границе 
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где о, - известное значение скорости, 


- на входной границе (границе втекания воздушного потока в помещение) 
Р=соп& (9) 
(условие Дирихле). 
Граничные условия для уравнения (3): 
начальные условия заключаются в задании распределения поля значений температуры в началь- 
ный момент времени (# = &), т. е. предшествующий расчетному. Они должны быть заданы в виде 
функции 


т =ЛСе,у,2) = (10) 
для пространственной задачи; 
- на боковых стенках 
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1 
данное уравнение обозначает, что со стенками помещения ведется теплообмен. Здесь должны 
быть заданы коэффициент теплоотдачи а и температура окружающей тело среды Т.; 
- на входной границе (границе втекания воздушного потока в помещение) 
Т =сопя- (12) 
известное и постоянное значение температуры втекаемого в помещение потока воздуха; 





к ОТ о 
- на выходной границе — =- т .(Т,), что означает, что поток температуры примеси за- 
т 


висит от её температуры, скорости о, и направления воздушного потока, индуцированного рабо- 
той вентиляции. 

В ходе экспериментальных работ граничные условия уточняются. 

Коэффициенты турбулентного обмена Ап. При турбулентном течении потока и турбулент- 
ном переносе коэффициенты Ап пропорциональны диссипируемой энергии и определяющему 
размеру рассматриваемого источника. Для вентилируемых помещений, зная величину энергии, 
поступающей в помещение, и определяющий размер источника, можно определить коэффициент 
турбулентного обмена в плоскости (х, У), м?/с: 


А=А,=0,25.58 1%, (13) 


где = - кинетическая энергия воздушных потоков в помещении; И = определяющий размер по- 


1 


мещения, равный для прямоугольных отверстий и сечений неправильной формы 
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и =УЕ ,м, (14) 
где Е— площадь сечения помещения, перпендикулярная направлению движения воздушных потоков. 
Коэффициент турбулентного обмена по высоте помещения находится из зависимости, 
предложенной М. Я. Берляндом для наибольших классов устойчивости воздушной среды: 


Ан = А’ `(2[2 1), (15) 
что допустимо для исследуемого помещения. Здесь А._, — коэффициент турбулентного обмена на 


высоте 1 м, 2, = 1м, 2- текущая высота помещения, м. 
Кинетическая энергия воздушных потоков в помещении = определяется из формулы 


&=8нс + 2 +8ди вс, (16) 
куда входят соответственно энергии приточных струй =‹, тепловых струй =, движущихся 
предметов =„ и вытяжных струй =,с . Величины &пс, &к, &ли Учитываются в модели рассеива- 


ния ВВ и находятся из формул: 
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где [/Г, - кратность воздухообмена в помещении свободным объёмом Г,; » - средняя ско- 
рость выхода воздуха из приточных отверстий, м/с; & — коэффициент местного сопротив- 
ления на выходе из воздухораспределителей; О’ - теплонапряжённость объема, Вт/м3; д = 


9,8 м/с?; 2— расстояние от полюса до рассматриваемой точки, м; п = 0,8 — эксперименталь- 
ная константа; С, — теплоемкость воздуха на удалении от источника, ; Т, — температура 


воздуха в помещении на удалении от источника; р, — плотность воздуха на удалении от ис- 
точника, & — коэффициент аэродинамического сопротивления движущегося предмета; Р’ — 
площадь поперечного сечения движущегося предмета; у, — скорость движущегося предме- 
та; х — средняя продолжительность движения предмета с данной скоростью в течение 1 ч; 
С, — масса воздуха в объеме помещения; 


так как С’, =И -р, где И- объем помещения, то: 
Е Ем [9 ь 
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Характеристика источника выделения вредностей 





В уравнение (1) входит О/(!) -— интенсивность выброса 1-го вредного вещества от /-го ис- 


точника. В нашем случае источник загрязнения -— это двигатель автомобиля ВАЗ 2104, который 
работает на холостом ходу. При этом с отработанными газами через шланговый отсос в воздух 
гаража поступают токсичные компоненты. Как показано в работе [5], одним из основных токсич- 
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ных компонентов, который учитывается при расчёте потребных воздухообменов на вентиляцию, 
является оксид углерода. 

Для определения его массового выброса необходимо экспериментально определить кон- 
центрацию СО в отработавших газах (ОГ) источника и измерить расход отработавших газов авто- 
мобиля. Для одиночного стационарного источника функция ОХ) может быть определена по фор- 
муле 


Сл .О 
Ск ОГерД 
© 3600 


где Сл — концентрация ВВ в ОГ от /-го источника загрязнения, г/м3; Оогер д - Расход отрабо- 


г/с, (21) 





тавших газов двигателя автомобиля на холостом ходу, м3/С. 
Явное теплопоступление в исследуемое помещение учитывается выражением 


О, =Сл, Оу, (То; -Т,)- Ал, кВт, (22) 
где Съог — изобарная теплоемкость ОГ на режиме холостого хода, Дж/кг ' К; Т р — темпера- 


Тура ОГ на холостом ходу, К; Т, — температура газовоздушной среды в помещении, К; 
А 


Д 
Экспериментальные исследования. С целью обоснования адекватности разработанной мате- 
матической модели распространения угарного газа, избытков теплоты и подвижности воздуха 
внутри заранее выбранного модельного объекта с простой геометрической формой и известным 
стационарным источником загрязнения были проведены экспериментальные исследования. 

Одновременно в заданные моменты времени замерялись концентрации угарного газа, 
температура газовоздушной смеси и скорость движения воздуха непосредственно в разных точках 
объёма исследуемой газовой камеры. Для оценки пространственного распределения вредных ве- 
ществ по помещению применялся стационарный метод отбора проб (функции устройств для от- 
бора и оперативного анализа проб выполняли электрохимические датчики газового анализа, ко- 
торые с прибором размещались в фиксированных точках помещения). 

Натурные исследования процессов проводились в гаражном помещении стандартных раз- 
меров (длина - 6,40 м, ширина - 3,30 м, высота - 2,30 м), схема которого (вид сверху) представ- 
лена на рис. 1. 
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Рис. 1. Схема размещения измерительных точек 
в гаражном помещении 


Измерения концентраций и температур газовоздушной среды за разные промежутки вре- 
мени проводились в пяти заданных точках камеры, где указаны их геометрические координаты: 
«а» (1.10,1.10,1.5), «Б» (1.10,2.20,1.5), «с» (3.20,1.65,1.5), «а» (5.20,1.10,1.5), «е» 
(5.20,2.20,1.5). Геометрические параметры источника загрязнения соответствуют координатам 
точки «с» — центр помещения (см. рис. 1). 

Согласно разработанной «Программе испытаний», замеры концентраций угарного газа в 
воздухе газовой камеры осуществлялись на высоте 1,5 м в характерных точках, представленных 
на схеме исследуемого объекта (см. рис. 1). 
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Измерения концентраций угарного газа в воздухе газовой камеры проводились в холод- 
ный период года (ноябрь) в течение часа по следующему графику (табл. 1). 






































Таблица 1 
Последовательность взятия проб в течение часа 

№ п/п Номер пробы на газовый анализ Время, прошедшее с начала эксперимента, с 
1 1 0 

2 П 175 

3 ш 325 

4 М 475 

5 \ 625 

6 МТ 1260 

7 УП 2100 

8 УШ 2700 

9 14 3300 

10 Хх 3600 











Для полноты исследования физической картины образования и распространения темпера- 
турных полей и полей концентраций ВВ в камере испытания проводились для трёх характерных 
случаев: 1) закрытые ворота (вентиляционные отверстия заклеены); 2) открыта одна створка во- 
рот; 3) ворота приоткрыты (щель размером 0,47 м). 


Приборное оснащение эксперимента. Результаты измерений 


Для измерения концентраций ВВ применялся газоанализатор Огадег Х-ат 5000 (Герма- 
ния): метод измерения — электрохимический; диапазон показаний — СО: 0-2000 ррт, №02: 
0-50 ррт, $502: 0-50 ррт (ррт — мг/м3); предел допускаемой основной погрешности - 5%. 

Микроклиматические условия в помещении (температура и подвижность воздуха) измеря- 
лись согласно схеме, представленной на рис. 1, на высоте 1,5 м. Для измерений скорости движе- 
ния воздушных потоков и температуры газовоздушной среды использовался высокоточный диф- 
ференциальный манометр Тез 521. Он имеет внутренний сенсор с измерительным диапазоном 
от 0,1 до 100 ГПа и является прибором для измерений скорости потока от 0,1 до 100 м/с. Приме- 
няется совместно с трубкой Пито. Модель Тезю 521-1 со встроенным сенсором давления имеет 
погрешность 0,2% от полной шкалы, а модель Тез 521-2 со встроенным сенсором давления 
имеет погрешность 0,1% от полной шкалы (рис. 2). 
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Рис. 2. Динамика изменения концентраций СО в газовой камере (точки а-е см. рис. 1) 
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Измерения температуры производились в холодный период года внутри помещения на высоте 
1,5 м. Использовались ртутные термометры с точностью измерений до 0,1 °С. Зависимости темпера- 
туры газовоздушной смеси в камере в заданных измерительных точках представлены на рис. 3. 
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Рис. 3. Динамика изменения температуры газовоздушной среды в камере (точки а-е см. рис. 1) 


Экспериментальное определение характеристик источника заг| рязнений 


Массовый выброс СО от источника загрязнений (.,(Г) определялся экспериментально. 


Для проведения измерений использовался 4-компонентный газоанализатор СО-СН-СО?-02 
“АВТОТЕСТ-01.03-М” с вычислением коэффициента избытка воздуха Л. Краткие технические ха- 
рактеристики анализатора газа АВТОТЕСТ 01.03 М: диапазон измерения содержания оксида угле- 
рода СО 0-7 об. %, абсолютная погрешность = 0,2%, относительная погрешность = 6%. 


Экспериментальное определение коэффициентов турбулентного обмена 


Как следует из выражения (13), определяющее влияние на коэффициенты турбулентного 
обмена кроме геометрических размеров помещения оказывает кинетическая энергия газовоздуш- 
ных потоков в помещении. Точное определение данных коэффициентов невозможно, так как в 
формулы (17) - (19) входят трудноопределимые экспериментально величины: кратность вентиля- 
ции помещения; коэффициент аэродинамического сопротивления и ряд других. Поэтому в задачи 
экспериментальных исследований входила возможность оценки вышеприведённых величин. 

Для корректного определения энергии приточных струй (17) в газовой камере были про- 
ведены экспериментальные работы по определению кратности вентиляции помещения гаража 
для двух исследуемых случаев: 1) открыта одна створка ворот; 2) ворота приоткрыты (щель раз- 
мером 0,47 м). С этой целью вначале создавалась максимальная загазованность в помещении, а 
затем контролировалось время “спада” концентрации СО до фоновых значений. Параллельно из- 
мерялась скорость вытекания газовоздушной смеси из камеры. Зная площади поперечных сече- 


ний отверстий вытекания газов „5 (2,56 м? и 1 м?), определяли потребный воздухообмен и крат- 
ность вентиляции (табл. 2). 
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Таблица 2 
Потребный воздухообмен и кратность вентиляции 
Исследуемый случай Время проветрива- Скорость вытекания Потребный воздухо- Кратность воздухооб- 
ния, Тр, Ч газов из камеры, ®„, | ви а м ть 
м/с (м/ч) ГИ у, 
Одна створка открыта 0,37 1 (3600) 3410 75 
Щель размером 0,47 м 0,52 1,5 (5400) 2700 58 

















Затем по формулам (14) - (19) были рассчитаны энергии приточных и тепловых струй с 
учётом определяющих размеров помещения для каждого из случаев (табл. 3). 


Энергии приточных и тепловых струй 


Таблица 3 





Исследуемый случай 


Энергия приточных 


струй, пс, м?/с 


Энергия тепловых 


струй, &тс , м’/с* 


Суммарная энергия в 
камере 5 , м?/с3 


Определяющий раз- 
мер помещения, 























1=УЕ,м 
Одна створка открыта 0,75 0,000227 0,750227 1,6 
Щель размером 0,47 м 0,58 0,000227 0,580227 1,0 





Наконец, в соответствии с формулами (4), (13) и (15) определили значения коэффициен- 
тов турбулентной диффузии и турбулентной теплопроводности в помещении для исследуемого 
угарного газа для трёх изучаемых случаев. Результаты расчётов сведены в табл. 4. 




















Таблица 4 
Коэффициенты турбулентной диффузии и теплопроводности 
Исследуемый Коэффициенты турбулентной | Коэффициенты турбулент- Коэффициент Коэффициент турбу- 
случай (вариант | диффузии помещения камеры, | ной теплопроводности по- турбулентной лентной теплопровод- 
расчёта) м?/с мещения камеры, Вт/м`К диффузии газа, ности газа, Вт/м"К 
М^/с 
Апх = Апу Ап Ях =Яу Л Осо Й,со 
Закрытые ворота 0 0 0 0 15,210 26,2:103 
Открыта одна 0,42 п 504 = 15,2108 26,2:103 
створка ворот Апг-1,5 $3 1,5 360 
= 0-3 
ром 0,2085 Е = 0,3 240 Яна = 360 15,210 26,2'10 























Коэффициенты турбулентной диффузии угарного газа Д., рассчитаны по формуле [6]: 
—7 3/2 
р- 4,3.107.ТГ 





— р.(у 


А 


| 1 1 
; + , 
+) М, М, 


(23) 


где М„, М», - молярные массы СО и воздуха соответственно, г/моль; у,, у, — мольные объёмы 


СО и воздуха, см3; р- давление газовоздушной смеси, Па; Т-— абсолютная температура га- 
зовоздушной смеси, К. 
Результаты модельных расчётов и сравнение их с экспериментальными данными. 
В связи с тем, что на практике в различных помещениях всегда работает общеобменная или ак- 
тивная вентиляция, то случай 1 (закрытые ворота — вентиляционные отверстия заклеены) являет- 
ся наиболее подходящим для теоритических исследований процессов диффузии газа, без учета 
вентиляционных потоков. Данный случай помогает выяснить влияние газов на диффузионную 
составляющую уравнений системы, описывающей процесс тепломассопереноса, поэтому было 
принято решение подробнее рассмотреть данный вариант исследования. 
Результаты модельных расчётов концентраций и температур газовоздушной среды внутри 
гаражного помещения для первого случая — закрытые ворота (вентиляционные отверстия заклее- 
ны) — представлены на рис. 4-7. 
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Рис. 4. Динамика изменения концентраций СО в газовой камере (точки а-е см. рис. 1) 
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Рис. 5. Динамика изменения температуры газовоздушной среды в камере (точки а-е см. рис. 1) 


Как следует из графиков зависимостей (см. рис. 4, 5), процесс роста температуры в замк- 


нутом объёме в значительной степени инерционный по сравнению с быстрым ростом концентра- 
ции угарного газа. 


Также следует отметить, что концентрации ВВ быстро выравниваются по величине в раз- 
ных частях исследуемого помещения (см. рис. 4), в отличие от температур газовоздушной смеси, 
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где имеет место значительное расхождение между параметрами источника (точка «с» — макси- 
мальная температура) и другими расчётными точками (см. рис. 5). Очевидно, что последнее об- 
стоятельство связано со сложными процессами теплообмена газовой камеры с окружающей средой. 

Картина сформированных полей концентраций ВВ и температур газовоздушной среды ис- 
следуемого помещения в вертикальном «срезе» на момент времени 2700 с представлена на 
рис. 6, 7. Как видно из результатов модельных расчётов, по высоте помещения имеет место чёт- 
кое расслоение концентраций СО, причём их максимум приходится в центре помещения у потол- 
ка. Тёплый угарный газ поднимается вверх, скапливаясь у потолка. В результате концентрации 
отличаются у пола и потолка в 2 раза. Несколько иное расслоение имеет место по температурам 
в помещении. Картина имеет симметричную форму, влияние на которую оказывает теплообмен с 
окружающей камеру средой. Необходимо отметить, что исследования проводились в «холодный» 
период года и разница температур ОГ и атмосферного воздуха составляла около 40-50 °С. 
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Рис. 6. Поле концентраций СО в газовой камере на высоте 1,5 м от пола 
(точки а-е см. рис. 1) 
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Рис. 7. Поле температур газовоздушной среды в камере на высоте 1,5 м от пола 
(точки а-е см. рис. 1) 
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Наконец, на рис. 8, 9 даны графики изменения концентраций и температур внутри закры- 
того помещения камеры во времени, полученные как экспериментально, так и в результате рас- 
чёта по разработанной математической модели. 
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Рис. 8. Сравнение теоретически рассчитанных концентраций СО 
в газовоздушной среде камеры по модели с экспериментальными данными 
(точки а, с см. рис. 1). Случай - ворота закрыты 
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Рис. 9. Сравнение теоретически рассчитанных температур 
газовоздушной среды камеры по модели с экспериментальными данными 
(точки а, с см. рис. 1). Случай — ворота закрыты 


Как видно из рис. 8, отклонение экспериментальных данных от расчётных концентраций 
СО в воздушной среде исследуемой камеры составляет не более 2-3%, что соответствует пределу 
допускаемой основной погрешности измерительной аппаратуры — 5% в интервале времени 0-175 с. 
Дальнейшее расхождение расчётных и экспериментальных зависимостей в интервале 175—625 с 
связано с тем, что газоанализатор Огадег Х-ат 5000 имеет ограниченный предел измерений по 
СО до 2000 ррт и не позволяет определить дальнейшее возрастание концентраций. 


872 


Вестник ДГТУ. 2011. Т. 11, № 6(57) 








На рис. 9 представлены графики динамики изменения температуры газовоздушной среды 
исследуемой камеры в течение часа, полученные как в результате численного моделирования, 
так и при проведении натурного эксперимента. Сравнение произведено для двух характерных 
точек помещения. 

Как следует из рис. 9, отклонение экспериментальных данных от расчётных температур 
газовоздушной среды камеры составляет не более 2-3 °С, что является вполне допустимым для 
моделирования таких достаточно сложных явлений тепломассопереноса веществ в замкнутом 
пространстве. 


Выводы. 

1. Разработана математическая модель процессов распространения газообразного веще- 
ства и избытков теплоты в воздушной среде гаражного помещения. 

2. Доказана адекватность представленной математической модели путём сравнения ре- 
зультатов численных расчётов искомых параметров внутри исследуемого помещения с парамет- 
рами, определёнными экспериментально на основе построенной физической модели. 

3. Показана возможность применения разработанной модели при исследовании загазо- 
ванности производственных помещений и для определения параметров микроклимата в рабочих 
зонах. 
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